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ABSTRACT 



We have searched for corotations (CR) in three Southern grand design spi- 
ral galaxies: NGC 1365, NGC 1566 and NGC2997. We have also introduced a 
method of quantifying errors in the phase diagram used to detect CR. 

We established the m=2 pattern CR at 12.1 kpc, 9.4kpc and 7kpc, for 
NGC 1365, NGC 1566 and NGC 2997, respectively. By using published rotation 
curves, we could determine spiral pattern angular speeds of 25.0km/sec/kpc, 
12.2km/sec/kpc and 17.6km/sec/kpc, respectively. 

A three armed component has been detected in NGC 2997, with the CR placed 
at 8.7 kpc with a pattern angular speed ^cRa = 12.7km/sec/kpc. 

An m=l component was detected in NGC 1566. We warily locate the CR 
at 7.1 kpc, with a pattern angular speed ^cr — 16.6km/sec/kpc. This pattern 
does not present ILR. 

Ages have been determined by studying the radial density profile of the m=2 
Fourier components in g (newly formed stars) and i (perturbing SDW supported 
by the disk of old stars), aided by the global aspect of the real spiral pattern 
in comparison with numerical simulations. The pattern is ~1200Myr old in 
NGC 1365, ~800Myr old in NGC 1566 and younger than SOMyr in NGC 2997. 

Subject headings: Galaxies: spiral, structure, kinematics and dynamics, 
individual (NGC 1365, NGC 1566, NGC 2997) - Methods: numerical 

1. Introduction 

Within the framework of the Density Wave Theory (Lin & Shu 1964), the spiral arms of 
disk galaxies are the manifestation of traveling waves. In a first approximation, the pattern 
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speed of a density wave is constant over the disk. It is widely known that galactic disks 
present differential rotation and in this way, in a given galaxy, there is a radial position 
where the density wave pattern speed coincides with the angular velocity of the stellar disk. 
This radial position is know as the corotation resonance (CR) radius. The density wave 
pattern is slower than the disk material at radii smaller than the CR, and it is faster at 
radii larger than the CR, inducing shocks at different arms side inward and outward the 
CR. Other resonances can be present in a given galaxy, mainly the ILR (Inner Lindblad 
Resonance), and the OLR (Outer Lindblad Resonance). These resonances are placed when 
the condition Qp = Q{R) ± k,/2 is satisfied. Here, Qp is the pattern speed, Q{R) is the 
angular velocity of the disk, k, is the epicyclic frequency, and the signals '+' and '— ' refer 
to the ILR and the OLR, respectively. 

Schweizer (1976) and Beckman & Cepa (1990) have previously discussed what would 
be the behavior of the colors across spiral arms if a shock generated by a spiral density 
wave (SDW) induces star formation. The main observable characteristic of this scenario 
are steeper azimuthal profiles and bluer color indices on the side where the shock front is 
located. 

Puerari & Dottori (1997) proposed a photometric method to detect CR based on the 
comparison of blue and infrared frames azimuthal profiles Fourier transform, under the 
basic idea that a shock-induced star formation in a SDW scenario produces an azimuthal 
gradient of ages across the spiral arms that has opposite signs on either side of the CR. 

The behavior of the azimuthal phase 0(r) was schematically exemplified in Fig. 2 of 
Puerari & Dottori (1997), for leading and trailing waves of S- and Z-type. Elmegreen, 
Elmegreen & Montenegro (1992) pointed out that evidences for CR are clearly seen in 
gas rich galaxies in the form of sharp end points to star formation ridges and dust lanes 
in two armed spirals. Photometric methods are widely used to detect resonances, mainly 
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due to the spar in observing time. Canzian (1998) rise some doubts on the abihty of 
photometric methods to uniquely determine CR in MIOO. For this galaxy, Elmegreen et 
al. (1989), using computing enhanced imagery, placed the CR radius at llOarcsec. This 
value was confirmed by Sempere et al. (1995), through the HI velocity field. On the other 
hand, two different methods based on Ha velocity field place the CR of M 100 at around 
75arcsec (Arsenault et al. 1988; Canzian & Allen 1997). That shows for MIOO that even 
the kinematical determination of CR from the velocity field of two different components do 
not lead to a unique result. Perhaps, the HI streaming motions over M 100 arms (Knapen 
et al. 1992) contribute to disturb SDW kinematics, indicating that perturbations in disks of 
galaxies sometimes produce more complex spatial structures and/or kinematical patterns 
than those suggested by simple models. 

In this paper, we use and discuss again the azimuthal profiles phase method (Puerari 
& Dottori 1997) to analyze the grand design spiral galaxies NGC 1365, NGC 1566 and 
NGC2997, recalling the curious result on multiple CR obtained by these authors for 
NGC 1832 and NGC 7479. NGC 1365 and NGC 1566 belong to Elmegreen & Elmegreen 
(1987) arm class 12 ('two long symmetric arms dominating the optical disk'), while 
NGC 2997 belongs to the arm class 9 ('two symmetric inner arms; multiple long and 
continuous outer arms'). 

Furthermore, in this paper we present a method to quantitatively determine errors. 
This procedure allows to determine a range of validity in radius for the azimuthal profiles 
phase method, a critical issue in Puerari & Dottori (1997) Letter. We apply the azimuthal 
profiles phase method on the rough g and i images, as done by Puerari & Dottori (1997), but 
also to images computer-treated with the method of Elmegreen, Elmegreen & Montenegro 
(1992), and that based on modes search by two-dimensional Fourier analysis (Kalnajs 1975; 
Considere & Athanassoula 1982; Puerari & Dottori 1992). 
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2. Data and Analysis 

2.1. Observations 

The galaxies discussed in this paper were observed in November, 1998, at the CTIO 
0.9m telescope. For each galaxy, three images with each of Gunn's g and i filters were 
taken, with exposure time 500 seconds each image. The standard reduction was carried out 
using the IRAF package^. Since our focus is morphology, no calibration were required. As 
a first step, all images were centered using field stars. Thereafter, field stars were removed 
using the IRAF imedit task. In order to enhance the perturbation over the galaxy's disks 
we subtract the mean radial light profile and then we normalize the rms variation of the 
intensity at each radius to a constant value throughout the image, as done by Elmegreen, 
Elmegreen & Montenegro (1992). The galaxies were then deprojected running a program 
that maps ellipses into circles. 

2.2. The Azimuthal Profile Phase Method 

Shock induced star formation in a stellar density wave scenario produces an azimuthal 
spread of ages across the spiral arms. At the CR, the angular velocity of the perturbation 
(Qp) and that of the stellar disk (Q) coincide. A comoving observer at the CR will see 
outward and inward the shock front to change from one side of the spiral to the other, 
consequently reversing the order in which young and older disk stellar population appear 



^The IRAF package is written and supported by the IRAF programming group at the Na- 
tional Optical Astronomy Observatories (NOAO) in Tucson, Arizona. NOAO is operated by 
the Association of Universities for Research in Astronomy (AURA), Inc. under cooperative 
agreement with the National Science Foundation (NSF). 
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in azimuthal profiles across the arms. In order to detect the shock front jump, Puerari & 
Dottori (1997) proposed to analyze the relative behavior of the SDW and shock front phases 
respectively, by means of the Fourier transform of azimuthal profiles lr{0) given for m=2 by 



F2{r) = I \r{0)e-^'^de 

J —-K 

The phase 0(r) can be obtained for m=2 as 



where Re and Im mean the real and imaginary parts of the complex Fourier coefficients. 
In Figs. 1 and 2 of Puerari & Dottori (1997), one can find a graphic representation of 
0(r) behavior for 'S' and 'Z' arms of trailing and leading character. This method was also 
successfully used by Aguerri, Beckman & Prieto (1998) to find CR in barred galaxies. 

The SDW phase and the newly formed stars one are respectively identified with 
0d«)2 = 02i and 0s/2 = 02^) where the right terms are obtained from infrared i and blue 
g images, respectively. An Azimuthal Profile Phase Difference (hereinafter, APPD) 
[©23 — 02i] = 0, will indicate the presence of a CR. 

We also search in this paper for the existence of the m=3 component. For this 
component the phase G is given by 

where 
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2.3. The Errors 

The basic idea for errors determination is that if the phase differences in the comparison 
of g and i images reflect an statistical fluctuation, it should also be reflected when different 
single filter images are compared among them. In order to quantitatively evaluate the errors 
intrinsic to the APPD method, we obtained three different images in g and i filters. We 
co-added gX^gl and g2+g3 and compute the phase difference of the resulting images. The 
same procedure is applied to the i images. When comparing single filter images of a given 
galaxy, the global aspect of the APPD is that of a high frequency statistical fluctuation 
(white noise), while APPD of g against i images produces low frequency phase differences, 
as we will see further on. We assume that an APPD that locally reach values larger than 
30" constitutes a real signal. 

2.4. Separation of m=2 and ni=3 Components 

In order to separate m=2 and m=3 spiral pattern components we apply two methods, 
namely, that of Elmegreen, Elmegreen & Montenegro (1992) and the one based on 
two-dimensional Fourier transform (Puerari & Dottori 1992). 

2.4- i- Elmegreen, Elmegreen & Montenegro Method (EEM92) 

These authors propose to separate the two- and three-fold symmetric part of a 
galaxy spiral pattern, S2 and S3 respectively, by making successive images rotations and 



subtractions. If /(r, 6) is the original image in polar coordinates, then 

S2{r,e) = I{r,e) - [I{r,e) - I{r,e + n)]T 
and 



Ssir,e) = 2Iir,e) - [/(r, ^) + I{r,e + ^)]T - [/(r, ^) - I{r,9-^)]T 



where T stands for truncation, meaning that pixels with negative intensities are set to 0. 



2.4-2. The Fourier method 

This method has been extensively discussed in a number of papers (Kalnajs 1975; 
Considere & Athanassoula 1982; lye et al. 1982; Puerari & Dottori 1992; Puerari 1993). In 
the Fourier method, an image is decomposed on a basis of logarithmic spiral of the form 
r = roexp( — — ^). The Fourier coefficients A{p,m) can be written as 

A{p,m) = — / / I{u,9)exp[—i{m9 + pu)]dud9 

^ J —TT J —OO 

Here, m = In r, r and 9 are the polar coordinates, m represents the number of the 
arms, p is related to the pitch angle a of the spiral as a = atan(— m/p), and I{u,9) is 
the distribution of light of a given deprojected galaxy, in a In r versus 9 plane. D is a 
normalization factor written as 



/ I{u, 9)dud9 

-n J —oo 
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Obs.: In practice, the integrals in -u = In r are calculated from a minimum radius (selected 
to exclude the bulge where there is no information on the arms) to a maximum radius 
(which extends to the outer limits of the arms in our images). 

The inverse Fourier transform can be written as 



Amd 



m 

where Pm stands for the radial density profile of the m component, written as 



Pm{n) = ^^P GUp)Aip,m)e'P- dp 



and Gm{p) is a high frequency filter used to smooth the A{p, m) spectra at the interval ends 
(Puerari & Dottori 1992), which has the form 



Gm{p) = exp 



I fp-pT '^ 



^max 



25 



where p^^x i^ ^^^ value of p for which the amplitude of the Fourier coefficients for a given 
m is maximum. The chosen interval ends (p+ = 50 and p_ = —50), as well as the step 
dp = 0.25, are suitable for the analysis of galactic spiral arms. 

So, by using 2D Fourier transforms, we can separate any m-mode. For example, the 
m=2 mode can be analyzed using S2{u,6) = S2{u)e'^^^ , and so on. 

3. Spiral Structure Properties in NGC 1365 

The disk subtracted and rectified to face-on images of NGC 1365 in the g and i filters, 
are presented in Figure 1. Two dimensional Fourier analysis shows that the spiral structure 
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in NGC 1365 is mainly composed of m=2 components, in both, g and i frames (Figure 2). 
The strong signal appearing in m=4 and m=6 are only aliases of m=2 (Puerari & Dottori 
1992). Nevertheless, the radial dependence of the m=2 arm intensity P2(r) is not similar in 
both frames, as Figure 7 shows. Up to lOkpc, the SDW intensity, represented by the old 
stellar component, P2i(r), predominates on that of the perturbed material, represented by 
the newly formed stars, P2g(r). From 13kpc outwards, P2j(^) weakens compared to 'P2g{r). 
When compared with numerical models, it becomes an interesting method to determine the 
pattern age, as we will see in section 3.2. 

3.1. CR Resonance in NGC 1365 

The identification of resonances in barred spirals is a matter of increasing interest for 
many reasons related to the SDW theory. By example, the self-consistent modeling of real 
systems Lindblad, Lindblad & Athanassoula (1996), the study of the time evolution of SDW 
(formation, dumping, etc) and related induced star formation (Junqueira & Combes 1996; 
Sellwood & Sparke 1988; Rautiainen & Salo 1999), and the study of non-linear coupling 
among different spirals modes (Tagger et al. 1987; Sygnet et al. 1988; Masset & Tagger 
1997b), are interesting, even open, problems. 

Resonances in NGC 1365 have been kinematically analyzed by Jorsater & van Moorsel 
(1995) and Lindblad, Lindblad & Athanassoula (1996). They place the CR resonance 
at 13.8 kpc and 15.2 kpc (Rcij/Rfoar 1-15 and 1.3), respectively. Lindblad, Lindblad & 
Athanassoula (1996), point that in models where the disk suffers a pure bar perturbation, 
a pattern speed higher than 20 ± 1 km/sec/kpc (corresponding to smaller CR radii) tends 
to give spirals more tighty wound than observed. They also point out that the CR radius 
should be even more external, if the perturbation is composed by a bar plus spiral arms. 
Nevertheless, as mentioned by these authors the evolution of the gaseous and/or stellar 
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system to its present form is beyond the scope of their paper. This evolution precisely 
transforms dramatically the spiral pattern parameters (number of arms, pitch angle, etc.), 
distinctly in old and young components, as can be seen in Junqueira & Combes (1996) 
self-consistent numerical experiments. This might explain the difference between Jorsater 
& van Moorsel (1995) and Lindblad, Lindblad & Athanassoula (1996) results. Our results 
are more similar to those of Jorsater & van Moorsel (1995), as we will see in the next 
paragraph. 

In Figure 3 we present the rectified g image, in a logr vs. 9 diagram (Elmegreen, 
Elmegreen & Montenegro 1992). In this plot we can see more precisely the different arm 
structures. The bar extends up to 10.3 kpc and its axis is strictly straight, as can be better 
seen in the western arm, which is less disturbed by the dust lane. The dust lane presents an 
open spiral form, with pitch angle a = 67° with extension similar to that of the bar. The 
bisymmetrical arms present two different pitch angles. At the beginning it has a = 16°, 
and extends itself from the end of the bar up to 17.5 kpc. The rest of the arm, up to 
^28 kpc, has a = 40°. It has to be pointed out that the pitch angle change is not an 
artifact introduced by the galaxy warping, since the galaxy was corrected by an inclination 
uj = 40° and the warping does not depart more than 15° from this value. Moreover, the 
warping begins at ?»25kpc (Jorsater & van Moorsel 1995). 

We present in Figures 4, 5 and 6 the behavior of the m=2 APPD of the original 
images ([02c, — Q2i]orig), and those transformed by EEM92 ([62c, — 02j]eem) and Fourier 
([©23 — 02i]FoMr) methods. The corresponding errors in g and i, are obtained according to 
the prescriptions of Section 2.3. 

In Figure 4, one can see that inside 5 kpc the S/N is too high, mainly in the i band, as 
to advance any conclusion. Between ^ 6 kpc and 15 kpc, [O23 — 02i] comfortably overcome 
the noise by many a. Beyond 15 kpc, it is again strongly affected by the noise. The CR 
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at the extreme of the bar is set at lO.Tkpc by this method. The more internal cuts are 
probably produced by the influence of the dust lane, whose curvature is different from that 
of the stellar bar. 

[Qsf — Qdw]eem and [Qsf — Qdwjpour prcscut a similar behavior than that of the 
original images, but with a better S/N. The CR is set at 13.8 and 11.8 kpc respectively 
by these two methods, with a mean value in close agreement to Jorsater & van Moorsel 
(1995) determination. Both methods show a tendency to unravel a CR at 25 kpc, distance 
at which the disk begin to warp. But the error in the i image increase strongly, and prevent 
any conclusion on the reality of this CR. Deeper imagery could help to better understand 
the SDW behavior in this region. 

According to the rotation curve of Persic & Salucci (1995), the mean value for the CR 
radius of rcj:j=12.1kpc leads to a pattern speed r2pattern=25.0Km/s/kpc. 

We have searched for the presence of m = 3 component, by means of EEM92 and 
Fourier transform methods, but this component is not significant in NGC1365 (see the 
spectrum for m=3 in Fig. 2). The m=l component appears to weak in the Fourier analysis, 
and prevent any CR study. 

3.2. NGC 1365 Spiral Pattern Age 

Junqueira & Combes (1996) performed numerical models of disk galaxies with n = l 
Toomre's disks of stars and gas. In all the experiments they used a radial scale length, 
a^ = 3.5 kpc and 4.0 kpc, for the star disk, and a^ = 6.0 kpc for the gaseous disk. The extend 
of their stellar and gaseous disks is 12 kpc and 16 kpc, respectively. Their models were 
completed with a Plummer's bulge with scale length, a^ ranging from 1/4.5 to 1/3 of as, and 
mass equal to that of the stellar disk. These models do not have halo. Junqueira & Combes 
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(1996) present the evolution of the stars and gas spiral patterns between 20Myr and 
2000 Myr for ten different experiments. NGC 1365 compares very well to their experiment 
E. It has a^ = 3 kpc, and a5 = l. In this numerical experiment, they found that the arm 
structure, which is multiple armed at 200 Myr, becomes two-armed at 800 Myr. From 
there on, the radial behavior of the m=2 Fourier intensity of young P2g and old P2J stellar 
population, which were fairly similar at 200 Myr, begin to differentiate between them. P2i 
predominates in the inner 1/4 of the disk length, and P2g outwards. We compare P2g and 
P2J in Figure 7, where we plotted corresponding Fourier radial intensities for perturbing and 
perturbed material (Junqueira & Combes 1996), after 1200 Myr (their Figure 10). As can be 
seen, the qualitative behavior is quite similar to that of NGC 1365, giving to this diagram 
a good chance as a model dependent spiral pattern age indicator. These experiments are 
specially suited for NGC 1365 because this galaxy does not possess an important halo, as 
clearly demonstrated by the highly Keplerian behavior of the rotation curve for r > Ri 20 kpc 
(Jorsater & van Moorsel 1995). 

4. Spiral Structure Properties of NGC 1566 

The images in the g and i filters for NGC 1566, disk subtracted and rectified to face-on, 
are presented in Figure 8. 

The rotation curve of this galaxy (Persic & Salucci 1995) presents a high dispersion, 
because it is projected practically face-on on the sky . For this reason, a photometric 
analysis of the resonances is very appropriated in the case of NGC 1566 (Elmegreen & 
Elmegreen 1990). 

The HII regions spiral pattern of this galaxy was already analyzed using Fourier 
transform by Puerari & Dottori (1990). The predominant pattern of this young stellar 
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population is m = 2, with m = 1 being progressively more important at radii larger than 
17kpc. The present analysis of g and i imagery shows up the same property (Figure 9). 
Contrary to the case of NGC 1365, in NGC 1566 the m = 2 SDW intensity in the old stellar 
component P2i does not predominates in any part of the disk on that of the newly formed 
stars, P23 (Figure 14). Both functions decrease similarly and fade up simultaneously at 
about 22kpc. 

For the one armed component, we note that between 10 and 13kpc Pij dominates on 
Pig, while in more external regions the contrary happens (Figure 15). 

4.1. CR Resonance in NGC 1566 

The resonances of NGC 1566 were studied by Elmegreen & Elmegreen (1990). The 
method they used was an identification of optical features together with constrains set by 
plausible rotation curves. They fitted radii for the ILR, OLR, CR, and other resonances. 

As in the case of NGC 1365, in Figures 11, 12, 13 we present the behavior of the APPD 
for the three methods used in this analysis. As before, the corresponding errors in g and i, 
are obtained according to the prescriptions of Section 2.3. 

Each one of the methods respectively places the CR at R~ 9.7 kpc (Fig. 11), Rsi 
8.5 kpc (Fig. 12), and Rr^ 10 kpc (Fig. 13) in good agreement among them. This distance 
corresponds to the CR assumed by Elmegreen & Elmegreen (1990). While some other 
intersections can be seen in the diagrams, the S/N indicates that their evidence is marginal. 

In Fig. 10 we present the rectified image in a logr vs. 9 diagram. As we can see, at the 
CR circle the main arms begin to thick outwards. The same behavior can be seen in the logr 
vs. 9 diagram of the symmetrized S2 image (not presented here). Other important property 
of this galaxy is that the arms, although ill defined between 10 and 15 kpc, are logarithmic 
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between 5 and 22kpc, with a pitch angle a ~ 36°. A short secondary arm, parallel to the 
main one, appears between 12 and 22 kpc. It is better defined in the SW arm, appearing 
as a rather asymmetric structure. Surprisingly a similar structure appears at 800 Myr in 
Junqueira and Combes experiment F2 (their Fig. 19), in agreement with the age determined 
in the next section for the SDW in this galaxy. The arms bend themselves between 20 and 
24 kpc, with a different pitch angle, although the feature is too weak to further analysis. 
We search also for m=3 component in NGC 1566, but it is not important. On the other 
hand, an m=l component is present, although tenuously defined. One can warily locate its 
CR at 7.1 kpc, without ILR, and pattern angular speed ^cr = 16.6km/sec/kpc. 

4.2. NGC 1566 Spiral Pattern Age 

As in the case of NGC 1365, we compared the spiral pattern properties with Junqueira 
& Combes (1996) models. Odd SDW components in Junqueira & Combes (1996) models 
appear at earlier evolutionary stages of the perturbation. Later on, the model spiral 
structure transforms itself into a two armed one. The more suitable model for the case of 
NGC 1566 is that of their experiment F, with n = 1 Toomre's disks of stars and gas, and 
radial scale length of the star and gaseous disks, ag = 4.0 kpc, and a^ = 6.0 kpc, respectively. 
The extend of the stellar and gaseous disks is 12 kpc and 16 kpc, respectively. The 
Plummer's bulge scale length for this model is a^ = as /4, and the bulge mass is equal to 
that of the stellar disk. 

The behavior of P23 and P2i (Figure 14) as well as the Fourier images of old and young 
stars are compatible with an SDW age of about 800 Myr in this galaxy. We do not know 
how important the halo is in this galaxy, in order to evaluate if the models are so suitable 
as in the former case. Nevertheless, Parkes 21 cm observations (Becker et al. 1988) indicate 
that the HI extends up to 90% of the optical disk, pointing also in this case to a not too 
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much extended halo. 

5. Spiral Structure Properties of NGC 2997 

We present in Figure 16 the disk subtracted and deprojected to face-on image of 
NGC2997. 

This galaxy is not as symmetric as NGC 1365 and NGC 1566. As discussed before, 
NGC 2997 is classified by Elmegreen & Elmegreen (1987) as arm class 9 ('two symmetric 
inner arms; multiple long and continuous outer arms'), while NGC 1365 and NGC 1566 are 
both arm class 12 ('two long symmetric arms dominating the optical disk'). 

Fourier analysis of HIT regions was performed by Considere & Athanassoula (1982), 
Puerari & Dottori (1992) and Garcia-Gomez & Athanassoula (1993). All studies show a 
predominance of m=2 and m=3 modes. At some radii, m=3 is even stronger than m=2. 
Near- infrared images of NGC 2997 (Block & Puerari 1999), reveal a strong asymmetry on 
the main inner arms, suggesting the presence of an m=l mode, besides the m=3 mode, in 
the older disk population (Block et al. 1994). The Fourier bisymmetric component in a K' 
image shows a pitch angle a = 25° . 

In Figure 17 we present the Fourier spectra for the g and i images of NGC 2997. The 
m=2 mode is, as expected, the predominant one. It is important to note that the m=l 
component, which appears at the noise level in the g image, presents a reasonable weight 
with respect to m=3 in i image, although weaker than m=3. Since in K', m=l dominates 
on m=3 (Block & Puerari 1999), and m=3 has practically a color independent weight 
relative to m=2, it is strongly suggesting that the m=l component is the only one really 
affected by dust obscuration. 
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5.1. CR Resonances in NGC 2997 

It was very difficult to apply the Phase method to the original images of NGC 2997. 
The m=2 phase did follow the internal arms, but becomes very chaotic at the distance 
of the arm bifurcation. The presence of several strong interarm features also contributes 
to increase the noise and makes the results less reliable. So, the methods that involve 
image processing, namely EEM92 and Fourier transform, are more suitable in the case of 
NGC 2997. 

To treat properly the main m=2 component, the part of the bisymmetric image in 
Figure 19 containing the secondary m=2 component in the E-W direction, from 7kpc 
outwards, was erased. Figure 20 shows the APPD for the main m=2 component. There is 
a clear cut at 7kpc, with very good S/N. The CR is located where the EW extension of the 
m=2 arms appears (Figure 19). This can be easily seen in Figure 18, where we plotted the 
symmetrized image in log r versus 9. From the rotation curve of Sperandio et al. (1995), we 
derived a 2-arms SDW angular velocity of ^sdw2 = 18km/s/kpc. The same place for the 
m=2 CR position is found from Fourier analysis (Figure 21). 

EEM92 method clearly reveals in this case the presence of a 3-arms component 
(Figure 23), each arm covering an arch of ^180° between 3kpc <r <22kpc. Similar 
3-arms morphology has been detected by EEM92 in NGC 628, NGC 1232, NGC 5457 and 
NGC 6912. The global sense of winding of this component is the same as that of the 2-armed 
one. Nevertheless, in the case of NGC 2997 the 3-arms are too wide at the beginning and 
present in this region a remarkable concavity of leading nature between 7kpc<r < llkpc. 
We dettached the concavity with a filled curve in Figure 23. The blow-up in the figure 
shows details of the ridge of one of the 3-arms concavity. APPD diagram (Figure 24), shows 
a intersection at 8.5 kpc, where the phase shift crosses downward, in agreement with the 
leading character of the described 3-arms concavity, indicating that star formation occurs at 
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the ridge crests. The detection of the elusive 3-armed SDW CR resonance in this galaxy is 
one important point of this study. Nevertheless, a better understanding would be necessary 
of what the phase diagram should do for a 3-armed component. 

From the rotation curve of Sperandio et al. (1995) we obtained a pattern speed for the 
3-arms SDW ^sdws ~ 13.7km/s/kpc. 

In Figure 18 we present NGC 2997 rectified image of the m=2 component as derived 
from EEM92 method, in a logr vs 6 diagram. The bifurcation of the m=2 component 
begins at the CR radius (rcr = 7kpc). The main m=2 spiral extends up to r=16kpc and 
the secondary between 7kpc and 21kpc. Both spirals are logarithmic, with pitch angle 
a = 15°. 7 and P=39°.4, respectively. 

5.2. NGC 2997 Spiral Pattern Age 

We compare the images and m=2 Fourier density profiles of NGC 2997 (see Figure 22) 
with Junqueira & Combes (1996) models. Here, the presence of asymmetries and higher m 
modes betrays the young age of the perturbation in NGC 2997. As discussed above, odd 
SDW components in Junqueira & Combes (1996) models appear at earlier evolutionary 
stages of the perturbation. In this case, it is difficult to assign an specific model from 
Junqueira & Combes (1996) which fits NGC 2997. Models E, F, G, and H show similar 
features at stages younger than 80 Myr. The m=2 Fourier radial density profiles present two 
differentiated parts: Inside 14kpc P23 ~ P2jj while outwards the old population function 
predominates on that of the young one, probably indicating that the external parts of the 
arms may have been produced by an older SDW. This last assessment has to be taken 
with caution, since Junqueira & Combes (1996) model F that fits this galaxy, seems to be 
oversimplified to explain the scenario of the SDW in NGC 2997. 
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6. Discussion and Conclusions 

The one-D Fourier transform method, apphed to g (blue) and i (infrared) images 
proposed by Puerari & Dottori (1997) to detect CR, has been applied to three grand design 
galaxies, namely, NGC 1365, NGC1566 and NGC2997. Some improvements have been 
introduced in this paper, to make the analysis more reliable. The original treatment on 
rough images of Puerari & Dottori (1997), was this time complemented with the analysis of 
computer treated images, with 2-D Fourier transform and EEM92 methods. This complete 
procedure presents complementary aspects, which gives advantages with respect to the 
original one. For example, to work on the original image is very good for pure m=2 arms. 
When other components appear, Fourier transform and EEM92 methods work better. 
Fourier method gives smaller errors than EEM92 method, but EEM92 method is more 
realistic than the former, since it is not constrained to describe the spiral as pure logarithmic 
ones. EEM92 method also have detected nicely the presence of the m=3 component in 
NGC 2997. This method is extremely efficient in detecting this component, while Fourier 
methods has to add many components. So, we confident can say that m=3 is not present 
in NGC 1365 and NGC 1566, and also that NGC 1566 presents an m=l component. 

The method presented here to determine errors in the APPD diagram, also represents a 
considerable improvement with respect to the original paper (Puerari & Dottori 1997), that 
strongly helps to identify which is the level of confidence in a CR determination. In turn, it 
also helps to establish the correct observational time to verify the dubious cases. From this 
point of view, it will be important to check Puerari & Dottori (1997) result on multiple CR 
in the galaxies NGC 1832 and NGC 7479. Nevertheless, more images of these galaxies are 
necessary to obtain a S/N in APPD diagram. The present form of this diagram, as a phase 
difference rather than as the superposition of the phases in each color, is more instructive 
than the former one, specially to visualize the S/N. 
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The determination of the arms extend has been made via logr vs. 9 diagram and via 
rotation curve, once know the CR radius. 

The arm complexity as well as the relative arm intensity for perturbing SDW P2i(r) 
and perturbed material P2g(i') in the Fourier diagram are related to the perturbations age. 
Although model dependent, a comparison with suitable models, allows to establish a scale 
of ages for the galaxies analyzed here. That leads to think that these three objects present 
different phases of a single phenomenon, in spite of the strong bar in NGC 1365. These 
conclusions have to be further checked with more elaborated models and for a larger sample 
of galaxies. 
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(a) (b) 



Fig. 1. — NGC 1365 disk subtracted and rectified to face-on images according to procedure 
described in section 2.1. (a) g image, (b) i image. 
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Fig. 2. — Two-dimensional Fourier spectra coefficients A{p,m) of NGC 1365, obtained as 
described in section 2.4.2. Left, g image, and right i image. The p variable is related to the 
pitch angle a hj p = ~-^^^- 
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Fig. 3. — NGC 1365 g band rectified image in log r versus 9 diagram. The horizontal line 
shows the CR loci at 13 kpc. According to Elmegreen, Elmegreen & Montenegro (1992) (Here- 
inafter EEM92). 
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Fig. 5. — APPD [©23 — 02i] for NGC 1365 obtained from images symmetrized according to 
EEM92 Method (solid line). Broken line shows the noise in g color and dotted line shows 
that in i color. See sections 2.4.1 and 2.3. 
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Fig. 6. — APPD [O23 — 02j] for NGC 1365 obtained from the Fourier ni=2 modes (solid line). 
Broken line shows g noise and dotted line shows i noise. See sections 2.4.2 and 2.3. 
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Fig. 7. — Comparison of Fourier radial density functions in blue P2g and infrared P2i (see 

section 2.4.2) for NGC 1365. Inset shows SDW and gas Fourier radial density functions for 
Junqueira & Combes (1996) model E after 1200 Myr. 
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(a) (b) 



Fig. 8. — NGC 1566 disk subtracted and rectified to face-on images according to procedure 
described in section 2.1. (a) g image, (b) i image. 
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Fig. 9. — Two-dimensional Fourier spectra coefficients A{p,m) of NGC 1566, obtained as 
described in section 2.4.2. Left, g image, and right i image. The p variable is related to the 
pitch angle a hj p = — -^^^^. 
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Fig. 10. — NGC 1566 g band rectified image in log r versus 9 diagram. The horizontal line 
shows the CR loci at 9.4 kpc. According to EEM92. 
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11. — APPD [©23 — 02i] for NGC 1566 obtained from original images (solid line). Broken 
shows g noise and dotted line shows i noise. 
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Fig. 12. — APPD [©2^ — 02j] for NGC 1566 obtained from images symmetrized according to 
EEM92 Method (solid line). Broken line shows the noise in g color and dotted line shows 
that in i color. See sections 2.4.1 and 2.3. 
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Fig. 13. — APPD [©23 — 02i] for NGC 1566 obtained from the Fourier ni=2 modes (solid 
line). Broken line shows g noise and dotted line shows i noise. See section 2.4.2. 
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Fig. 14. — Comparison of Fourier radial density functions in blue P2g and infrared P2i (see 

section 2.4.2) for NGC 1566. 
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Fig. 15. — Comparison of Fourier radial density functions in blue Pig{r) and infrared Pii{r) 

(see section 2.4.2) for NGC 1566. 



-39- 



(a) (b) 



Fig. 16. — NGC2997 disk subtracted and rectified to face-on images according to procedure 
described in section 2.1. (a) g image, (b) i image. 
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Fig. 17. — Two-dimensional Fourier spectra coefficients A{p, m) of NGC 2997, obtained as 
described in section 2.4.2. Left, g image, and right i image. The p variable is related to the 
pitch angle ahj p = — -^^^. 
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Fig. 18. — NGC2997 g band symmetrized EEM92 S2 and rectified image in log R versus 9 
diagram. The horizontal line shows the CR loci at 7kpc. According to EEM92. 
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Fig. 19.— NGC2997 S2 symmetrized image according to EEM92 Method. 
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Fig. 20. — APPD [O23 — 02j] for NGC2997 obtained from images symmetrized according 
to EEM92 Method, inner arms between 4 < r < 11 kpc (solid line). Broken line shows the 
noise in g color and dotted line shows that in i color. 
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Fig. 21.— APPD [023 - 02i] for NGC2997 obtained from the Fourier m=2 modes (solid 
line). Broken line shows g noise and dotted line shows i noise. See section 2.4.2. 
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Fig. 22. — Comparison of Fourier radial density functions in blue P2g and infrared P2i (see 



section 2.4.2) for NGC2997. 
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Fig. 23. — NGC2997 S3 symmetrized image according to EEM92 Method. Curves drawn 
at the 3-arm beginning points to the leading fronts of star formation. The blow-up shows 
details of one of the three fronts. See section 5.1 paragraph 3. 



-47- 



~l 1 T" 




_l I l_ 



o 
I 5 10 15 

RADIUS (Kpc) 

Fig. 24. — APPD [©33 - Qsi] for NGC2997 obtained from m=3 Fourier modes (solid line). 

This Fourier analysis was performed in S3 symmetrized images according to EEM92 Method. 

Broken line shows the noise in g color and dotted line shows that in i color. 
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Table 1. Galaxy parameters. 



Galaxy Type RA Dec i° PA° v (km/s) Conv(kpc/") 



NGC1365 SBb(s)b 03h33m36.4s -36d08m25s 40 220 
NGC1566 SAB(rs)bc 04h20m00.6s -54d56ml7s 37 60 
NGC2997 SA(s)c 09h45m38.6s -31dllm25s 40 80 



1636 


0.1074 


1496 


0.0967 


1087 


0.073 
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Table 2. Resonances mean value. h=75 Km/s/Mpc 

Galaxy CR ILR OLR m QcR{km/ s/kpc) 

NGC1365^ 12.1±0.9 0.6 18.3 2 25.0±1.9 

NGC1566'' 9.4±0.4 3.2 14.6 2 12.2±0.9 

7.r •■• 14.8 1 16.6 

NGC2997<^ 7.0 2.7 9.7 2 17.6 

8.7 6.0 10.9 3 12.7 



'^The rotation curve was taken from Persic & Salucci (1995) 
and Jorsater & van Moorsel (1995). 

*^The rotation curve was taken from Persic & Salucci (1995). 
^One-arm CR was estimated from Log r vs O diagram of 

the antisymmetric images g and i according to Elmegreen, 

Elmegreen & Montenegro method. The S/N is low and the 

value should be taken with ware. 

'^The rotation curve was taken from Sperandio et al. (1995). 
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